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Executive Summary 

The objective of this initiative, funded by NASA’s Aviation Safety Program, is to model, validate, 
and predict, with high fidelity, the microstructural evolution of third-generation high-refractory Ni-based 
disc superalloys during heat treating and service conditions. This initiative is a natural extension of the 
D ARP A- AIM (Accelerated Insertion of Materials) initiative with GE/Pratt- Whitney and with other 
process simulation tools. Strong collaboration with the NASA Glenn Research Center (GRC) is a key 
component of this initiative and the focus of this program is on industrially relevant disk alloys and heat 
treatment processes identified by GRC. Employing QuesTek’s Computational Materials Dynamics 
technology and PrecipiCalc precipitation simulator, physics-based models are being used to achieve high 
predictive accuracy and precision. Combining these models with experimental data and probabilistic 
analysis, “virtual alloy design” can be performed. The predicted microstructures can be optimized to 
promote desirable features and concurrently eliminate nondesirable phases that can limit the reliability 
and durability of the alloys. The well-calibrated and well-integrated software tools that are being applied 
under the proposed program will help gas turbine disk alloy manufacturers, processing facilities, and 
NASA, to efficiently and effectively improve the performance of current and future disk materials. 
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Introduction 


In support of higher fidelity application of the DARPA-AIM acceleration methodology previously 
demonstrated on the IN100 and R88DT disc alloys, a collaborative model and database validation study is 
supported by the NASA Aviation Safety Program to better quantify uncertainty and improve prediction 
accuracy of the QuesTek PrecipiCalc (ref. 5) multiphase/multicomponent precipitation simulation code. 
The nominal compositions of relevant disk alloys are tabulated in Table 1. 


TABLE 1. — NOMINAL COMPOSITIONS (WT%) OF THE DISK ALLOYS STUDIED UNDER THIS PROGRAM 


at.% 

Ni 

Cr 

Co 

Mo 

W 

A1 

Ti 

Nb 

Ta 

C 

B 

Zr 

ME3 

Bal. 

13.1 

20 

3.8 

1.9 

3.5 

3.6 

1.1 

2.3 

0.04 

0.03 

0.05 

LSHR 

Bal. 

13 

21 

2.7 

4.3 

3.5 

3.5 

1.5 

1.6 

0.03 

0.03 

0.05 

Alloy 10 

Bal. 

10.2 

14.9 

2.7 

6.2 

3.7 

3.9 

1.9 

0.9 

0.03 

0.03 

0.1 

RR1000 

Bal. 

15 

18.5 

5 

___ 

3 

3.6 

___ 

2 

0.027 

0.015 

0.06 


Program Status 

The program is on schedule for this reporting period. Major accomplishments for this reporting 
period are detailed in the Gantt chart in Figure 1. 



Yr 1 

Yr 2 

Yr 3 

Q1 | Q2 | Q3 | Q4 

Q1 | Q2 | Q3 | Q4 

G1 | G2 | Q3 | G4 

Task 1 - Fundamental Multicomponent Models 

Literature Surrey 

Evaluation and Assessment of Thermodynamics/Mobility 
Molar Volume Model 




Task 2 - Development of Calibration and Validation Data 

Identification of Data and Design of Focused Experiments 
3D APT Experiment and Analysis 




Task 3 - y' Precipitation Modeling 

Determination of Modeling Approaches 
PrecipiCalc Calibration and Validation 


IS 


Spatial Prediction of Microstructure Evolution 



4 

Task 4 - Precipitation Modeling for Embrittling Phases 

Determination of Modeling Approaches 
Model Calibration and Validation 




Prediction of Embrittling Phase Formation 



4 

Task 5 - Software Implementation and Dissemination 

Software Development 

4 

► 


Software Dissemination 



4 

Task 6 - Microstructure Variation and Optimization 

Evaluation ofVariation Sources 

Microstructure Variation Analysis 

Strategies for Improving Life and Durability of Alloys 



4 

Task 7 - Project Management and Repotting 

4 

► 4 

► 4 


Figure 1 . — Project Gantt chart of the proposed work. 
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Task 1 — Fundamental Multicomponent Models 


Thermodynamics 

Our assessment indicates that Ni-DATA version 7 in combination with NIST-Mobility is most 
favorable, because it reproduces the y’ -phase precipitation behavior accurately in Ni-Al-Cr and Ni-Al-Cr- 
X (X = Re, W) model alloys, as well as in commercial disc alloys. Thermodynamic databases compared 
in this study include ThermoTech Ni-DATA versions 4 to 7 (ref. 12), Ni-NIST (not including Nb) 

(ref. 1), Thermo-Calc TCNI1 (not including Mo, Nb, and Ta) (ref. 3), and CompuTherm PanNickel 
(ref. 1 1) in combination with NIST-Mobility (ref. 2) and MOBNil (ref. 10). Compared to previous 
versions, Ni-DATA version 7 contains updates to improve the y’ -phase solvus temperature prediction for 
recently developed Ni-base alloys with high refractory elements. The equilibrium concentrations 
calculated with Ni-DATA version 7 are in good agreement with the equilibrium concentrations at infinite 
time. See Table 2 for the comparison summary for the Ni-Al-Cr system, where Ni7 was found to be best, 
according to the root mean square (RMS) analysis. 


TABLE 2.— COMPARISON AND ANALYSIS OF COMPOSITIONS FROM APT MEASUREMENT 
AND THERMODYNAMICS EQUILIBRIUM PREDICTIONS WITH VARIOUS DATABASES 

FOR A Ni-Al-Cr SYSTEM 



600 °C 


at.% 

Ni Al Cr 

APB 

RMS 

600 °C 
misfit 

y’ vol. frac. 

y 

Equilibrium matrix 

APT 
2 a 

81.26 

0.18 

3.13 

0.08 

15.61 

0.18 


0.14 



Ni5 

Ni6 

Ni7 

TCNI1 

Ni-NIST 

81.40 

81.40 

81.34 

81.17 

80.60 

3.64 

3.64 

3.62 

3.93 

5.20 

14.96 

14.96 

15.04 

14.90 

14.20 


0.82 

0.82 

0.75 

1.07 

2.50 



Y ? 

Equilibrium 
composition 
of y'-precipitates 

APT 

2a 

76.53 

0.50 

16.69 

0.44 

6.77 

0.30 

0.19 

0.01 

0.38 

-0.72% 

14.8% 

0.5% 

Ni5 

Ni6 

Ni7 

TCNI1 

Ni-NIST 

74.98 

74.98 

74.97 

75.59 

75.96 

16.23 

16.23 

16.21 

14.55 

12.25 

8.80 

8.80 

8.82 

9.57 

11.79 

0.19 

0.19 

0.19 

0.18 

0.15 

2.08 

2.08 

2.11 

3.52 

6.70 

-0.58% 

-0.58% 

-0.57% 

-0.75% 

-0.89% 

12.4% 

12.4% 

12.4% 

12.3% 

0.0% 

Critical y'-nucleus 
composition 

APT 

2a 

72.40 

2.20 

18.30 

1.80 

9.30 

1.40 

0.16 

0.03 

1.61 

-0.22% 


Ni5 

Ni6 

Ni7 

TCNI1 

Ni-NIST 

74.91 

74.91 

74.91 

75.36 

No 

17.37 

17.37 

17.40 

16.43 

precipital 

7.71 

7.71 

7.69 

8.21 

tion 

0.20 

0.20 

0.20 

0.19 

1.84 

1.84 

1.84 

2.16 

-0.52% 

-0.52% 

-0.51% 

-0.61% 



The comparisons in Ni-10Al-8.5Cr-2X (X = Re, W) model alloys again support Ni-Data version 7 as 
most accurate. In the Ni-Al-Cr-Re alloy, the phase relations measured by Yoon et al. (ref. 8) suggest a 
partial solute-trapping behavior of slow-diffusing species Re, see Table 3. 
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TABLE 3.— COMPARISON AND ANALYSIS OF COMPOSITIONS FROM APT MEASUREMENT 
AND THERMODYNAMICS EQUILIBRIUM PREDICTIONS WITH VARIOUS DATABASES 

FOR A Ni-Al-Cr-RE SYSTEM 



u 

O 

o 

o 

OO 


Ni 

Cr 

Al 

Re 

APB 

RMS 

misfit 

y’ vol. 











frac. 

y 

Equilibrium 

APT 

81.07 

10.04 

6.74 

2.15 






matrix 

2 a 

0.16 

0.04 

0.08 

0.02 


0.09 





Ni 7 

80.55 

8.94 

8.35 

2.16 


1.95 





Ni-NIST 

80.22 

8.79 

8.87 

2.12 


2.47 



Y’ 

Equilibrium 

APT 

76.17 

4.97 

18.05 

0.81 

0.19 


-0.71% 

25.4% 


composition 
of y'-precipitates 

2 a 

0.26 

0.10 

0.18 

0.04 

0.00 

0.21 


0.2% 


Ni 7 

75.29 

5.68 

18.80 

0.23 

0.20 

1.18 

- 0 . 68 % 

15 . 9 % 



Ni-NIST 

75.82 

6.13 

17.65 

0.39 


1.30 

-0.77% 

7.2% 


80% equilibrium 

APT 

76.33 

5.46 

16.92 

1.29 

0.17 



20.3% 


composition of 
y'-precipitates 

2 a 

0.30 

0.16 

0.26 

0.08 

0.00 

0.32 


0.3% 


Ni 7 

75.22 

5.56 

18.99 

0.23 

0.20 

2.33 


12 . 8 % 


TABLE 4.— COMPARISON AND ANALYSIS OF COMPOSITIONS FROM APT MEASUREMENT 
AND THERMODYNAMICS EQUILIBRIUM PREDICTIONS WITH VARIOUS DATABASES 

FOR ANi-Al-Cr-W SYSTEM 












y’ vol. 


800 °C 


Ni 

Al 

Cr 

W 

APB 

RMS 

misfit 

frac. 



APT 

81.31 

5.83 

11.52 

1.34 





y 

Equilibrium 

2 a 

0.07 

0.04 

0.05 

0.02 


0.07 



matrix 

Ni 7 

81.44 

6.22 

10.53 

1.80 


1.16 





Ni-NIST 

80.48 

8.57 

8.84 

2.11 


3.91 





APT 

76.30 

17.00 

3.91 

2.80 

0.21 


-0.49% 

38.2% 

Y' 

Equilibrium 

composition 

2 a 

0.08 

0.07 

0.04 

0.03 

0.00 

0.09 


0.1% 

Ni 7 

75.39 

18.00 

4.20 

2.42 

0.22 

1.11 

- 0 . 36 % 

32 . 1 % 

of y'-precipitates 



Ni-NIST 

75.87 

16.61 

6.42 

1.09 


3.06 

-0.61% 

18.7% 


Composition trajectories of isothermal y’ precipitation in simple Ni-Al-Cr-X alloys were examined. In 
a Ni-5.2Al-14.2Cr at.% alloy during isothermal precipitation at 1 1 12 °F (600 °C), Sudbrack et al. 

(ref. 7) at the Northwestern University Center for Atom Probe Tomography (NUCAPT) reported the y/y’ 
interfacial free energy as 22 to 23 mJ/m 2 . They also reported that the effective solute diffusivities are one 
order of magnitude smaller than the predicted intrinsic diffusivities in the matrix phase. Employing these 
values, the simulation results for the evolution of phase fraction and particle size are consistent with their 
Atom-Probe Tomography (APT) analysis. 

The time evolution of particle size measured by APT suggests the precipitation composition trajectory 
can be modeled as an unstable equilibrium with capillarity. The composition trajectory calculated by 
PrecipiCalc in the parallel tangent composition mass balance mode is consistent with the unstable 
equilibrium Thermo-Calc (ref. 13) calculations by adding capillary energy. However, our simulation 
results underestimate the Cr-enrichment in the y’ critical nucleus. As Cr reduces the interphase lattice 
misfit, an elastic energy could be added to the y’ free energy in our thermodynamic calculations to 
account for this. 

To further examine the solute-trapping behavior of slow-diffusing refractory elements such as Mo, 

Hf, Ta, W, and Re, we conducted a literature review on the diffusivity of these species. Karunaratne and 
Reed (ref. 4) reported interdiffusion coefficients estimated by diffusion profile measurements. In dilute 
additions, the interdiffusion coefficients display little concentration dependence, but the values can differ 
by an order of magnitude. A systematic trend is observed: elements furthest from the center of the 
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periodic table display the largest values and elements at the center the smallest. W and Re are most potent 
for retarding the y’ coarsening kinetics, and in turn most prone to kinetic trapping. 

However, Sudbrack et al. (ref. 6) report no such solute-trapping in Ni-10Al-8.5Cr-2W aged at 
1472 °F (800 °C) up to 264 hr. The prolonged aging seems to eliminate the trapping behavior. Also, we 
confirmed that prolonged aging eliminates trapping for W in ME3. Thus, there is no apparent need to 
account for solute-trapping in our thermodynamic calculations. 

The equilibrium concentrations in Ni-10Al-8.5Cr-2W calculated with Ni-DATA version 7 are in 
good agreement with the extrapolated equilibrium concentrations. Comparisons of y’ phase fraction and 
the at-temperature interphase lattice-parameter misfit with the database predictions for these alloys again 
support the Ni-DATA version 7 database as most accurate, and suggest the corresponding error in APB 
energy and misfit is acceptable within current structure/property model uncertainty. 

Mobility 

To assess accuracy of available atomic mobility databases for combination with thermodynamic 
databases to predict multicomponent diffusivities, linear diffusion multiples were prepared by Dr. Ivan 
Locci at the NASA Glenn Research Center. A 0.2-in. thick disk of pure Ni was diffusion bonded by hot 
pressing to a 0.125-in. thick disk of ME3 on one side and to a 0.125-in. thick of Alloy 10 on the other, for 
4 hr at 1610 °F in vacuum. After the initial bond was formed, the diffusion couple was subjected to 
additional annealing in a horizontal tube furnace in an argon atmosphere at two temperatures, 2000 °F 
(1093 °C) or 1700 °F (927 °C) for 100 or 300 hr. 

QuesTek’s DICTRA (ref. 9) simulations with Ni-DATA version 7 and NIST-Mobility of the 
diffusion multiple at the hot-pressing temperature 1400 °F predict no significant diffusion, as shown 
below in figure 2. Therefore, subsequent simulations are ignoring the hot-pressing step. 



Figure 2. — DICTRA predicted composition trajectories of the 
AlloylO-Ni diffusion couple after 1400 °F hot press. 
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Figure 3. — Comparison of microanalysis and DICTRA predicted composition traces of the ME3-Ni diffusion 
couples at different aging conditions. 
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Alloy-10 Mi AlloylO Mi 




Figure 4. — Comparison of microanalysis and DICTRA predicted composition traces of the AlloylO-Ni diffusion 
couples at different aging conditions. 


The NIST mobility database employed in the AIM project was compared with the Thermo-Calc 
MOB1 database in combination with various thermodynamic databases for ME3, Figure 3, and Alloy 10, 
Figure 4. In order to compare microprobe analyses of the diffusion multiple with DICTRA simulations, 
the average Matano interface (defined as the interface across which equal number of atoms have crossed 
in both directions) was equated to the origin of the calculated profiles. Assuming y’ precipitates act only 
as sink or source of solute for diffusion, i.e. no diffusion through y’, the agreement is reasonable for Co, 
Cr, Mo, and W in ME3/Ni diffusion couple held at 2000 °F (1093 °C) for 100 hr. Before setting the 
Matano interface as the origin, a significant discrepancy was noticed for W; however this was resolved 
with the Matano interface adjustment. With the exception of Al, the data does not confirm the predicted 
strong nonmonotonic profiles near the interface. 

Furthermore, the current comparison shows that the diffusion of Nb, Ta, Al, and Ti is less than 
predicted. Of these elements, Ta seems to suffer from microprobe detection efficiency. Nb, Al, and Ti are 
y’ -stabilizers and could be affected by the no-diffusion-through-y’ assumption. Dr. Misra from QuesTek 
Innovations will present these results at the upcoming NIST diffusion workshop. The results were then 
forwarded to Dr. Campbell at NIST for a mobility database calibration. If needed, the diffusion multiple 
will be heat-treated at a supersolvus temperature to eliminate complications from precipitates. 
Nonetheless, the predicted transition width between y+y’+MC and y is in excellent agreement with 
experimental measurement. 
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Molar Volume 


Comparison of relative interphase misfit provides a useful metric of database accuracy for the 
complex disc alloys. The molar volume model constructed in QuesTek’s prior AIM program was applied 
to calculate the lattice misfit in disk alloys. A model was implemented at QuesTek to calculate the 
nonmonotonic temperature dependency of lattice misfit, see Figure 5. The disc alloys fall into 2 groups of 
relatively high and low misfit. ME3, LSHR, and Alloy 10 show large positive misfit, consistent with the 
aligned precipitates observed in SEM. Of the third-generation disc alloys, only the RR1000 alloy is 
predicted to maintain the low misfit of the IN100 and R88DT alloys. In support of this prediction, while 
the fine intragranular precipitates in Alloy 10 show the cuboidal morphology and ordered arrays consistent 
with the predicted higher misfit, the y’ in RR1000 shows a spheroidal morphology and disordered 
distribution consistent with its predicted low misfit, see Figure 6. 



Figure 5. — Predicted yVy misfit at high and ambient temperatures 
based on y’/y compositions calculated by equilibrium. 
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Task 2 — Development of Calibration and Validation Data 

High Temperature Equilibrium Age Results 

To assess equilibrium phase fractions and compositions in commercial disc alloys, samples of ME3, 
LSHR, AlloylO and RR1000 disc alloys were given a supersolvus treatment at 2190 °F (1199 °C) 
followed by water quench, and subsequent 20~1000-hr treatments at temperatures of 2000 °F (1093 °C), 
1700 °F (930 °C), and 1400 °F (760 °C). Focal Electrode Atom Probe (FEAP) at Northwestern 
University was used to analyze the composition and fine y’ microstructure of selected samples. Fike a 
traditional 3-Dimensional Atom Probe (3DAP), the FEAP microscope by Imago Scientific Instruments 
achieves true 3 -dimensional atomic-scale analysis by using a high electric field to remove individual 
atoms from material surfaces and a position-sensitive detector to record information that reveals the 
atom’s position and identity. The incorporation of a local electrode eliminates or mitigates many of the 
performance limitations of traditional 3DAPs. The FEAP microscope at Northwestern University 
analyzes significantly larger volumes in much less time, compared to the traditional 3DAP. 

The first sample analyzed under FEAP was ME3 water quenched from 2190 °F (1 199 °C). In total, 
8.7 million atoms were analyzed with the new larger-FOV FEAP detector. The overall bulk composition 
is very close to the expected composition. The precipitates formed during quenching have a sphere- 
equivalent diameter of about 15 to 30 nm. They are enriched predominantly in Ni-Ti-Al-Ta, and also in 
Nb, B, C, and O. They are depleted predominantly in Cr-Co-Mo, and also somewhat in W, Fe, and Si. 
The measured compositions show good agreement with our multiphase PrecipiCalc simulation. 

In addition, the high-temperature matrix composition of ME3 aged at 2000 °F for 1000 hr was also 
analyzed by FEAP. 6.2 million atoms were analyzed with the new larger-FOV FEAP detector. The 
precipitates formed during quenching are enriched predominantly in Ni-Ti-Al-Ta, and also in Nb, B, C, 
and O. They are depleted predominantly in Cr-Co-Mo, see Figure 7, similar to the results from ME3 
water quench sample. Also, the slow-diffusing W showed solute-trapping behavior in y’ nanoscale 
particles. The measured compositions are compared with our multiphase PrecipiCalc simulation, see 
Table 5. 
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Figure 7. — y and y composition analysis of LEAP experiment on ME3 
2000 °F/1000 hr isothermal hold sample. Only fine y particles (formed 
during quench from 2000 °F) were observed in LEAP. 
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An experimental y’ and MC phase fraction at 2000 °F equilibrium were determined to be 24.8±1.5 
and 0.33 percent, respectively, by mass-balancing the LEAP measurement and calculated y’ and MC 
phase compositions. At bottom of Table 5, it is clearly seen that Ni-DATA version 7 is most favorable. 
Although the low temperature y’ equilibrium fraction are similar between several databases, the predicted 
y’ fraction at 2000 °F can vary by 7 percent (see Figure 8), though the predicted y’ compositions do not 
show large departure between databases (see Table 6). 


TABLE 5.— ME3 2000 °F MATRIX COMPOSITION COMPARISONS (INCLUDING THE FINE y’ FORMED 
DURING QUENCH) BETWEEN LEAP AND DATABASE PREDICTION 


y matrix 

at.% 


Ni 

A1 

Cr 

Co 

Ti 

Mo 

W 

Nb 

Ta 

C 

RMS 


Ni7 

46.4% 

6 . 0 % 

18.5% 

21.8% 

2 . 7 % 

3.0% 

0.7% 

0.4% 

0.5% 

0.01% 

2 . 18 % 


Ni5 

45.7% 

5.6% 

19.7% 

22.1% 

2.4% 

3.1% 

0.7% 

0.3% 

0.3% 

0.01% 

2.89% 


PanNickel 

44.3% 

5.3% 

20.3% 

23.0% 

2.4% 

2.9% 

0.8% 

0.4% 

0.5% 

0 . 05 % 

3.61% 


NIST-Ni 

44.7% 

5.1% 

19.9% 

23.8% 

2.4% 

3.2% 

0.8% 


0.2% 


3.81% 


LEAP 

44.9% 

6 . 7 % 

17.5% 

21.7% 

3 . 5 % 

3.4% 

1.0% 

0.5% 

0.5% 

0 . 04 % 



2a 

0.04% 

0.02% 

0.03% 

0.03% 

0.01% 

0.01% 

0.01% 

0.01% 

0.01% 

0.002% 


phase 
fraction at 
2000 °F 

Ni7 y' 

60.2% 

12.0% 

2.8% 

13.0% 

8.8% 

0.3% 

0.3% 

1.3% 

1.2% 

0.0% 


24.8+1.5% 

Ni7 MC 

0.0% 

0.0% 

0.2% 

0.0% 

24.5% 

0.2% 

0.2% 

10.5% 

14.5% 

47.9% 


0.33% 

System total 

49.6% 

7.5% 

14.5% 

19.5% 

4.3% 

2.3% 

0.6% 

0.7% 

0.7% 

0.2% 




Ni7 

NIST 

Ni4~6 

Pan- 

Nickel 

TCNI1 








y' phase 
fraction 

25.5% 

29.4% 

30.7% 

32.3% 

8.2% 









c 

o 

o 

po 


CD 

CO 

03 

sz 

cl 

03 

E 

E 

03 

CD 



Figure 8. — Calculated equilibrium y’ phase fraction versus temperature with different 
databases for ME3, showing large variation in 2000 °F prediction. 
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TABLE 6.— PREDICTED ME3 y’ COMPOSITIONS AT 2000 °F 


Database 

ME3 y’ site-fraction at 2000 °F (1093 °C) 

Ni-DATA 
version 7 

(Nio.8oCoo.i7Cr 0 .o2)3 (Al 0 . 48 Ti 0 . 35 Nb 0 . 05 Ta 0 . 05 Cr 0 . 04 Ni 0 . 02 ) 

Ni-DATA 
version 4~6 

(N i 0 .79Co 0 . 1 8 ^ 0 . 02)3 ( AI0.47T L34T a 0 . 06 Nb 0 . 06 Cr 0 . 04 Ni 0 . 02 ) 

PanNickel 

(Nio.82Co 0 .i7)3 (Al 0 . 49 Ti 0 . 33 Cr 0 . 06 Nb 0 . 05 Ta 0 . 04 Ni 0 . 0 iMO 0 . 0 i) 

Ni-NIST 

(Nio.86Co 0 .i2Cr 0 .oi)3 (Al 0 . 52 Ti 0 . 36 Ta 0 . 08 Cr 0 . 02 Ni 0 . 02 ) 


Using the total sample composition as representing the y composition at 2000 °F (1093 °C), the root 
mean square (RMS) between the calculated matrix composition vectors and the solute contents of the 
LEAP composition vector indicates that Ni-DATA version 7 gives the best agreement (see Table 5). The 
discrepancy is mainly caused by A1 and Ti, which are both underestimated in the calculations. It is 
noteworthy that the C composition prediction of PanNickel is closest to the measurement, which would 
result in a higher-fidelity description of the MC precipitation. 

The systematic underestimation of Ti in the calculated y matrix is also noted in QuesTek’s earlier 
simulations in DARPA-AIM initiative on commercial Ni-base superalloys, as shown in Table 7, and may 
have to be addressed in this program. 


TABLE 7.— COMPOSITION COMPARISONS FOR INI 00 AND RENE88DT STUDIED UNDER AIM PROGRAM 


at.% 

INI 00 (3 step) 

Ni 

A1 

Cr 

Co 

Ti 

Mo 

W 

Nb 

V 

Zr 

B 

C 

y Matrix 

PpC 

37.5 

2.7 

28.4 

26.1 

0.1 

3.8 



1.3 





3 DAP 

32.0 

2.1 

30.5 

30.8 

0.2 

3.4 



0.9 

0.007 

0.059 

0.011 

Primary y' 

PpC 

61.1 

14.1 

3.1 

12.1 

8.8 

0.4 



0.3 





EDS 

63.1 

12.4 

3.7 

11.1 

8.0 

0.8 



0.7 




Secondary y' 

PpC 

63.2 

15.3 

2.7 

10.7 

7.1 

0.5 



0.6 





3 DAP 

62.8 

14.7 

2.8 

10.1 

7.9 

0.9 



0.7 

0.003 

0.087 

0.012 

Tertiary y' 

PpC 

64.6 

17.2 

3.0 

9.3 

3.7 

1.1 



1.1 





3 DAP 

63.1 

16.3 

2.5 

9.1 

5.3 

1.9 



1.7 

0.052 

0.054 

0.000 

System total 

45.7 

10.0 

13.0 

17.0 

4.6 

1.9 



0.8 




Rene88DT 













(continuous cooling) 













y Matrix 

PpC 

50.1 

1.4 

26.5 

16.3 

0.3 

3.6 

1.6 

0.02 






3 DAP 

48.0 

1.3 

28.0 

17.0 

0.7 

3.4 

1.5 

0.01 



0.078 

0.046 

Secondary y' 

PpC 

68.0 

10.4 

2.0 

6.5 

11.2 

0.2 

0.6 

1.19 






3 DAP 

71.4 

10.3 

1.6 

4.2 

11.2 

0.5 

0.7 

0.01 



0.049 

0.004 

Tertiary y' 

PpC 

68.9 

11.2 

1.8 

5.8 

10.4 

0.2 

0.7 

1.03 






3 DAP 

69.4 

10.9 

2.9 

5.1 

9.3 

1.0 

1.3 

0.01 



0.064 

0.026 

System total 

51.6 

2.2 

24.6 

15.5 

1.1 

3.4 

1.6 

0.10 






To prepare the experimental results for PrecipiCalc calibration, the evolution of mean precipitate 
radius was assessed in Alloy 10, LSHR, and ME3. Eugine Kang, a student at the Northwestern 
University, conducted image analysis at GRC to measure the y’ equivalent-circle diameter from SEM 
images, and the results are summarized in Figure 9. The image analysis indicated a twofold increase of 
equivalent-circle diameter from intragranular to grain-boundary. In addition, the estimated ME3 y’ size 
after solution treat and water quench is about 20 nm in diameter, which is in good agreement with LEAP 
results of 15 to 30 nm discussed earlier. 
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Figure 9. — Summary of isothermal hold experimental 
results for ME3, LSHR and AlloylO. 


Single Sensor DTA Results 

For an effective experimental calibration of nucleation related material parameters, Dr. Alexandrov at 
the Ohio State University conducted Single Sensor Differential Thermal Analysis (SSDTA) 
measurements to detect the transformation onset temperatures in rapidly cooled arc-melted samples of 
ME3, LSHR, AlloylO, and RR1000. 40 samples of 1/8-in. diameter pins were analyzed. QuesTek 
processed the thermograms to reliably derive the y’ precipitation onset. 

The first set of SSDTA samples was complicated by the possible solidification segregation as samples 
were melted and solidified in the arc-melter. Therefore, we alternatively tested the second set using 
SSDTA in a furnace with Ar or He cooling to avoid melting. The results were still in good agreement 
with the button-melting experiments. The y’ precipitation onset was measured at 2053 °F (1123 °C) for 
the first population and at 1567 °F (853 °C) for the second population. The interfacial energy of 
intragranular y’ precipitates is being calibrated to these SSDTA furnace experiments. The results analyzed 
by the DTA analysis software are summarized in Table 8. A typical SSDTA curve is presented in Figure 
10 showing the first (A) nucleation burst. The transformation onset temperature is the “A s ” column, which 
will be used to calibrate y Vy surface energy in next section, while other identified temperatures are not 
currently used. 
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TABLE 8.— RESULTS SUMMARY OF OSU’S SSDTA FURNACE COOL SAMPLES 


[Temperatures are in °C.] 


Sample 

Cooling 

media 

CR1 160/1 125 

C/s 

As 

A F 

Additional 

effect 

Api 

B s 

B f 

A 10-2 

Ar 1 1 cf/h 

35 

1145 

1130 


989 

854 

750 

A10-3 

He 50 cf/h 

59 

1147 

1127 

1105-1069 




A 10-4 

He 50 cf/h 

63 

1145 

1127 

1121-1106 

964 

852 

712 

LSHR 1 

Ar 1 1 cf/h 

49 

1123 

1103 


961 



LSHR2 

Ar 1 1 cf/h 

39 

1123 

1103 


950 

853 

758 

LSHR 3 

He 50 cf/h 

68 

1126 

1110 

1091-1083 

950 



LSHR 4 

He 50 cf/h 

74 

1126 

1103 

-1091 

946 

880 

751 

ME3-1 

Ar 1 1 cf/h 

41 

1128 

1108 


893 

618 


ME3-2 

Ar 1 1 cf/h 

48 

1125 

1104 


899 

617 


ME3-3 

He 50 cf/h 

69 

1126 

1110 

1080-1050 

865 

587 


ME3-4 

He 50 cf/h 

70 

1126 

1103 

1087-1027 




RR1 000-1 

Ar 1 1 cf/h 

28 

1108 

1090 

1007- 

901 

805 


RR1 000-2 

Ar 1 1 cf/h 

35 

1108 

1095 

1013- 

928 

815 


RR1 000-3 

He 50 cf/h 

67 

1110 

1085 


914 

862 


RR1 000-4 

He 50 cf/h 

64 

1106 

1087 

1006- 

961 

884 



Data 


1160 


1100 


1040 


RR 12 Ar.red 

1160/1100: 2.14s, 27.97C/S 
RR 22 Ar.red 

1160/1100: 1.73s,34.78C/s 
RR 3-1 He. red 
1160/1100: 0.90s, 66.53C/S 
RR 4-1 He. red 
1160/1100: 0.93s, 64.46C/S 


0 1 
Time(secs) 


Figure 10. — Typical SSDTA curve for RR1000 showing 
primary burst event. 


Using LEAP, we also analyzed ME3 Ar gas-quenched from 2150 °F (1175 °C). This specimen came 
from the Single Sensor Differential Thermal Analysis (SSDTA) conducted at Ohio State University, 
which rapidly acquires the cooling profile. In total, 32 million atoms were analyzed. The overall bulk 
composition is very close to the expected composition. Larger y ? precipitates, about 40 to 60 nm in 
diameter, were observed together with smaller gamma’ precipitates, about 5 nm or less in diameter. Both 
sets of precipitates are enhanced predominantly in Ni-Ti-Al, and depleted in Cr-Co. The smaller 
precipitates showed a flat W composition profile, while the larger precipitates showed a distinct W 
depletion. A high-resolution SEM analysis on ME3 showed consistent results with the LEAP analysis, 
see, Figure 11. 
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Comparison of tomographic atom probe and high resolution SEM analyses of ME3 superalloy: 

Two distributions of y precipitates 




Figure 1 1 . — Result summary of LEAP and high resolution SEM experiments of ME3-1 Ar SSDTA sample. 

Task 3 — -y’ Precipitation Modeling 
Determination of Modeling Approaches 

The data developed in Task 2 are used in y’ precipitation modeling discussed in this section. A logical 
procedure of calibrating PrecipiCalc model parameters is summarized in the following table. Once the 
fundamental CALPHAD thermodynamics and mobility databases are calibrated through equilibrium age 
and diffusion couple experiment, SSDTA results of precipitation onset temperatures are used to determine 
the coherent surface energy (a CO h) and coherent elastic energy (G e i), which will be verified with non- 
isothermal quench results. Then the results of intragrain y’ from equilibrium age experiments will be used 
to calibrate incoherent surface energy (ai nco h) and interfacial dissipation parameters, while the results of 
grain boundary y’ will be used to calibrate diffusivity scaling factor (D sca i e /in C oh)* The coherent to 
incoherent transition will be justified with available literature models. 
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Experiments 

CALPHAD 

Thermodynamics and 
Mobility 

Precipitation Kinetics 
Model Parameters 

Equilibrium Age 
+ LEAP and EDS 

Matrix/Precipitation 
Compositions and Phase 
Fractions 


Diffusion Couple 
+ Microanalysis 

Dscale/coh 


SSDTA 


a coh 


G el 

Nonisothermal Quench 
+ LEAP 
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Equilibrium Age 
+ SEM/TEM for y’ size and 
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a incoh„ G scale/incoh 
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Figure 12. — Impact of diffusivity scaling factor to ME3-1 Ar 
SSDTA PrecipiCalc simulation. 

Calibration of PrecipiCalc With SSDTA Results 

SSDTA determines the onset of transformation, which is dominated by the nucleation. According to 
the nucleation theory, the most important material parameters affecting the nucleation rate are driving 
force (which is determined by the bulk thermodynamics), surface energy and elastic coherency energy — 
i.e., parameters affecting nucleation energy barrier. In this section, we present the results of calibrating the 
surface energy, while ignoring the coherency energy. 

Unlike the conventional approach using precipitation coarsening results to determine surface energy, 
the use of onset of transformation has a benefit to remove the dependency of diffusivity. As seen in 
Figure 12, varying the diffusivity scaling factor in PrecipiCalc does not affect the onset temperature, 
though it affects the final y’ particle sizes. Hence, our approach of calibrating surface energy with 
transformation onset temperature removes any uncertainty of diffusivity. 

The transformation onset temperatures from the SSDTA experiments were determined from a SSDTA 
data processing software. It takes the measured temperature versus time, fits the baseline prior to the 
transformation with a choice of function suitable for describing cooling, and then determines the 
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departure point of measured data from fitted function. To determine the onset temperature from 
PrecipiCalc simulation with physically equivalent interpretation as closely as possible to SSDTA, here 
are the steps (see Figure 13 as an example with ME3-lAr) 

(1) Process SSDTA measured temperature profiles to remove the latent heat contribution; 

(2) Calculate material compositions relevant to the y’ precipitation from matrix y phase — the high 
temperature phases (borides, carbides and undissolved y’) are removed with equilibrium calculations at 
highest SSDTA measured temperatures using Ni-DATA 7 databases. The undissolved y’ only occurs to 
Alloy 10 samples, and more discussion on this will follow. 

(3) Perform PrecipiCalc simulations with above information using Ni-DATA 7 and NIST mobility 
databases, with estimated surface energy; 

(4) Collect time (or temperature) evolution results of y’ volume fraction (see Figure 13(a)), and 
compositions of matrix and y’, compute the time (or temperature) evolution of molar enthalpy 

(see Figure 13(b)); 

(5) Compute temperature derivatives of the molar enthalpy with respect to temperature, dH/dT, (see 
Figure 13(c)), 

(6) Determine the transformation onset temperature where dH/dT changes by more than 
10 percent of average dH/dT values, with decreasing temperature. 



Figure 13. — Determination of transformation onset temperature 
(vertical red line) from PrecipiCalc simulation with ME3-1Ar 
using 0.0275 J/m 2 as an example. 
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Table 9 summarizes the PrecipiCalc results and comparison to the SSDTA measured transformation 
onset temperatures (2nd column), for all 4 alloy samples. We varied the surface energy from 0.0325 to 
0.0225 J/m 2 with 0.0025 J/m 2 decrement. The root mean square (RMS) of the difference between 
prediction and measurement is calculated for each material and estimated surface energy. 

TABLE 9.— CALIBRATION OF y’/y SURFACE ENERGY WITH THE TRANSFORMATION ONSET 


TEMPERATURE DETERMINED BY SSDTA 



SSDTA 
Onset T 
(°C) 

PrecipiCalc Simulation, Estimated Surface Energy (J/m 2 ) 

0.0325 

0.03 

0.0275 

0.025 

0.0225 

Onset 

T(°C) 

AT(°C) 

Onset 

T(°C) 

AT(°C) 

Onset 

T(°C) 

AT(°C) 

Onset 

T(°C) 

AT(°C) 

Onset 

T(°C) 

AT(°C) 

ME3-lAr 

1128 



1122.5 

5.5 

1127 

1 

1131.2 

-3.2 



ME3-2Ar 

1125 



1123.8 

1.2 

1126.8 

-1.8 

1130.9 

-5.9 



ME3-3He 

1126 



1121 

5 

1125.4 

0.6 

1130 

-4 



ME3-4He 

1126 



1119.6 

6.4 

1123.8 

2.2 

1127 

-1 



RMS 





4.94 


1.54 


3.94 



RRlOOO-lAr 

1108 



1096.6 

11.4 

1100.8 

7.2 

1105 

3 

1109.2 

-1.2 

RR1 000-2 Ar 

1108 



1096.9 

11.1 

1100.7 

7.3 

1104.5 

3.5 

1108.2 

-0.2 

RR1000-3He 

1110 



1094.5 

15.5 

1098.7 

11.3 

1103 

7 

1107.3 

2.7 

RR1000-4He 

1106 



1094.4 

11.6 

1098.6 

7.4 

1103.4 

2.6 

1106.7 

-0.7 

RMS 





12.53 


8.48 


4.39 


1.52 

LSHR-lAr 

1123 

1123.3 

-0.3 

1128.6 

-5.6 

1132.9 

-9.9 

1137.2 

-14.2 



LSHR-2Ar 

1123 

1123.6 

-0.6 

1128.9 

-5.9 

1133.6 

-10.6 

1137.9 

-14.9 



LSHR-3He 

1126 

1122.1 

3.9 

1127 

-1 

1131.9 

-5.9 

1135.9 

-9.9 



LSHR-4He 

1126 

1121.3 

4.7 

1126.3 

-0.3 

1131.3 

-5.3 

1135.4 

-9.4 



RMS 



3.07 


4.10 


8.27 


12.35 



A10-2Ar 

1145 



1136 

9.05 

1140.8 

4.2 

1144.7 

0.3 



A10-3He 

1147 



1139.1 

7.9 

1144.8 

2.2 

1149.1 

-2.1 



A10-4He 

1145 



1134.9 

10.1 

1140 

5 

1144.1 

0.9 



RMS 





9.06 


3.98 


1.33 




Based on this analysis, the optimal surface energies (shown with green background) can be 
determined with minimum RMS for each alloys. That is, 0.0275 J/m 2 for ME3, 0.0225 J/m 2 for RR1000, 
0.031 J/m 2 for LSHR (average of 0.0325 and 0.03), and 0.025 J/m 2 for Alloy 10. Note that minimum 
RMS for LSHR is higher than other alloys, because the experimental transformation onset temperatures 
for Ar and He cooling are in the opposite trend with conventional wisdom. Furthermore, it is observed 
that every decrement of 0.0025 J/m 2 in surface energy introduces an increase of about 4 -5 °C in RMS, 
which is about the accuracy of SSDTA thermocouple reported by OSU. Hence, the determined surface 
energy could have a ±0.0025 J/m 2 uncertainty, and then the surface energies of these four alloys are not 
that different from each other. However, the higher estimated surface energy of Alloy 10 could be due to 
the coherency elastic energy effect, which will be investigated in the 2nd year of this program. 

The concern on AlloylO SSDTA samples was that the unsolutionized large y’ could affect the 
subsequent nucleation of finer y\ The amount of undissolved y’ in AlloylO-2Ar SSDTA sample was 
estimated to be 9.5 percent fraction (see Figure 14), which is comparable to the predicted 8.4 percent 
equilibrium y’ fraction based on Ni7 database. The calculated onset temperatures of AlloylO presented in 
Table 9 were based on the assumption that undissolved y’ does not grow and does not affect the subsequent 
y’ formation, allowing us to remove the large y’ completely in those PrecipiCalc simulations. To confirm our 
assumption, we performed a Alloy 10-2Ar simulation with the presence of large undissolved y’ (with a mean 
size of about 3 pm and phase fraction of 8.4 percent). Following the quench in SSDTA, the large y’ grows 
slightly and reduces the available supersaturation. The transformation onset temperature was lower than 
without the large y’, but only by 1.6 °C (see Figure 15 showing the volume fraction evolution of the fine y’). 
Hence, we concluded the assumption to ignore large undissolved y’ is reasonable. 
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Figure 14. — NASA SEM observation of Alloy10-2Ar SSDTA 
sample, y’ volume fraction was estimated to be around 
9.5 percent. 



1060 1080 1100 1120 1140 1160 1180 

Temperature (C) 


Figure 15. — Comparison of the fine y’ fraction from two PrecipiCalc 
Alloy10-2Ar simulations with and without primary y\ showing the 
impact to transformation onset temperature is minimum. 


Calibration of PrecipiCalc With Isothermal Hold Experiments 

The microstructure of disc alloys aged at 2000 °F (1093 °C) for 1000 hr demonstrated the importance 
of treating both intragranular and grain boundary precipitation. PrecipiCalc contains both homogeneous 
and heterogeneous (dislocation and grain boundary) nucleation models. The simulations are being refined 
for separate nucleation and interfacial mobility parameters for grain boundary behavior in parallel with 
intragranular precipitation. Parameters specific to heterogeneous nucleation on grain boundary are grain 
diameter, wetting angle of precipitate interface on grain boundary, and external scaling factor to the 
nucleation rate to account for uncertainty. 

The work in this area is still in progress and will be reported in next report. 
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Preliminary Conclusions and Next Steps 

Results so far indicate the ThermoTech Ni-DATA version 7 thermodynamic database gives the best 
agreement with measured phase compositions in equilibrated microstructures, providing acceptable 
accuracy of predicted parameters such as APB energy and interphase misfit. More accurate predictions 
will likely require incorporation of coherency elastic energy for the high misfit third-generation disc 
alloys. The fact that the discrepancy of observed critical nucleus composition for the Ni-Al-Cr alloy is in 
the direction of lower misfit suggests the addition of elastic energy may improve prediction of y’ 
composition at the nanoscale. 

The technical project milestone in Task 1 of the first project year shown in Figure 1 has been 
achieved. The software development milestone in Task 5 (Figure 1) has also been achieved through the 
demonstration of suitable PrecipiCalc simulation for available data. 

Further PrecipiCalc simulations will employ separate nucleation and interfacial mobility parameters 
for grain boundary behavior in parallel with intragranular precipitation. 

Future validation studies will include carbide evolution and long-term stability with respect to 
undesirable TCP phases. 

Dr. Olson from QuesTek is planning to attend Superalloys 2008 and prepared an extended abstract 
titled “Precipitation Model Validation in 3rd Generation Aeroturbine Disc Alloys.” 

References 

1 . Database provided by Dr. Ursula Kattner at NIST. 

2. Campbell, C.E., W.J. Boettinger, and U.R. Kattner, 2002. Acta Materialia 50: 775-792. 

3. TCNI1 — TCS Ni-based Superalloys database, http://www.thermocalc.com/Products/ 

Databases/TCNI 1 .htm. 

4. Karunaratne, M.S.A. and R.C. Reed, 2003. Interdifusion of the platinum-group metals in nickel at 
elevated temperature. Acta Materialia , 51:2905-2919. 

5. Jou, H.-J., P. Voorheese and G.B. Olson, 2004. Computer simulations for the prediction of 
microstructure/property variation in aeroturbine disks, Superalloys 2004 , p. 877. 

6. Sudbrack, C.K., D. Isheim, R.D. Noebe, N.S. Jacobson, and D.N. Seidman. 2004. The Influence of 
Tungsten on the Chemical Composition of a Temporally Evolving Nanostructure of a Model Ni-Al- 
Cr Superalloy. Microsc. Micro anal. 10: 355-365. 

7. Sudbrack, C.K., R.D. Noebe, and D.N. Seidman. 2007. Compositional pathways and capillary effects 
during isothermal precipitation in a nondilute Ni-Al-Cr alloy. Acta Materialia. 55: 119-130. 

8. Yoon, K.E., R.D. Noebe, and D.N. Seidman. 2007. Effects of rhenium addition on the temporal 
evolution of the nanostructure and chemistry of a model Ni-Cr-Al superalloy. II: Analysis of the 
coarsening behavior. Acta Materialia. 55: 1 145-1157. 

9. DICTRA, Diffusion Controlled TRAnsformation simulation tool http://www.thermocalc.com/ 
Products/Dictra.html 

10. TCS Ni-alloys mobility database 

http://www.thermocalc.com/Products/Databases/Descriptions/DBD_MOBNil.pdf. 

11. http://www.computherm.com/databases.html, PanNickel database from CompuTherm LLC, Madison, 
WI. 

12. Saunders, N. in “Superalloys 1996” eds. R. Kissinger et al. TMS, Warrendale PA, 1996, p. 115. 
(http://www.thermotech.co.Uk/databases.html#Ni-DATA) 

13. Thermo-Calc software tool, http://www.thermocalc.com/Products/TCC.html. 


NAS A/CR— 2008-2 15199 


19 



REPORT DOCUMENTATION PAGE 

Form Approved 
OMB No. 0704-0188 

The public reporting burden for this collection of information is estimated to average 1 hour per response, including the time for reviewing instructions, searching existing data sources, gathering and maintaining the 
data needed, and completing and reviewing the collection of information. Send comments regarding this burden estimate or any other aspect of this collection of information, including suggestions for reducing this 
burden, to Department of Defense, Washington Headquarters Services, Directorate for Information Operations and Reports (0704-0188), 1215 Jefferson Davis Highway, Suite 1204, Arlington, VA 22202-4302. 
Respondents should be aware that notwithstanding any other provision of law, no person shall be subject to any penalty for failing to comply with a collection of information if it does not display a currently valid OMB 
control number. 

PLEASE DO NOT RETURN YOUR FORM TO THE ABOVE ADDRESS. 

1. REPORT DATE (DD-MM-YYYY) 
01-04-2008 

2. REPORT TYPE 

Final Contractor Report 

3. DATES COVERED (From - To) 

4. TITLE AND SUBTITLE 

Microstructure Modeling of 3rd Generation Disk Alloys 


5a. CONTRACT NUMBER 

NNC07CB01C 






5b. GRANT NUMBER 






5c. PROGRAM ELEMENT NUMBER 

6. AUTHOR(S) 

Jou, Herng-Jeng 




5d. PROJECT NUMBER 






5e. TASK NUMBER 






5f. WORK UNIT NUMBER 

WBS 698259.02.07.03.03.02 

7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 

QuesTek Innovations LLC 
1820 Ridge Avenue 
Evanston, Illinois 6020 1 

8. PERFORMING ORGANIZATION 
REPORT NUMBER 

E- 16471 

9. SPONSORING/MONITORING AGENCY NAME(S) AND ADDRESS(ES) 

National Aeronautics and Space Administration 
Washington, DC 20546-0001 


10. SPONSORING/MONITORS 
ACRONYM(S) 

NASA 






11. SPONSORING/MONITORING 
REPORT NUMBER 

NASA/ CR-2008-2 15199 

12. DISTRIBUTION/AVAILABILITY STATEMENT 

Unclassified-Unlimited 
Subject Category: 26 

Available electronically at http://gltrs.grc.nasa.gov 

This publication is available from the NASA Center for AeroSpace Information, 301-621-0390 

13. SUPPLEMENTARY NOTES 

14. ABSTRACT 

The objective of this initiative, funded by NASA's Aviation Safety Program, is to model, validate, and predict, with high fidelity, the 
micro structural evolution of third-generation high-refractory Ni-based disc superalloys during heat treating and service conditions. This 
initiative is a natural extension of the DARPA-AIM (Accelerated Insertion of Materials) initiative with GE/Pratt- Whitney and with other 
process simulation tools. Strong collaboration with the NASA Glenn Research Center (GRC) is a key component of this initiative and the 
focus of this program is on industrially relevant disk alloys and heat treatment processes identified by GRC. Employing QuesTek’ s 
Computational Materials Dynamics technology and PrecipiCalc precipitation simulator, physics-based models are being used to achieve 
high predictive accuracy and precision. Combining these models with experimental data and probabilistic analysis, “virtual alloy design” 
can be performed. The predicted microstructures can be optimized to promote desirable features and concurrently eliminate nondesirable 
phases that can limit the reliability and durability of the alloys. The well-calibrated and well-integrated software tools that are being applied 
under the proposed program will help gas turbine disk alloy manufacturers, processing facilities, and NASA, to efficiently and effectively 
improve the performance of current and future disk materials. 
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